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Abstract 



We have detected new HD absorption systems at high redshifts, z a b s = 2.626 and z a b s = 1.777, identified in 
the spectra of the quasars J 0812+3208 and Q 1331+170, respectively. Each of these systems consists of two 
subsystems. The HD column densities have been determined: logNn D = 15.70 ± 0.07 for za = 2.626443(2) and 
logNg D = 12.98 ± 0.22 for z B = 2.626276(2) in the spectrum of J 0812+3208 and logNg D = 14.83 ± 0.15 for 
z c = l. 77637(2) and logNg D = 14.61 ± 0.20 for z D = l. 77670(3) in the spectrum of Q 1331+170. The measured 
HD/H2 ratio for three of these subsystems has been found to be considerably higher than its values typical 
of clouds in our Galaxy. We discuss the problem of determining the primordial deuterium abundance, which 
is most sensitive to the baryon density of the Universe fli,. Using a well-known model for the chemistry of a 
molecular cloud, we have estimated the isotopic ratio D/H=HD/2H2 = (2.97±0.55) x 10~ 5 and the corresponding 
baryon density SI,, ft 2 = 0.0205 l^ao - This value is in good agreement with Q b h 2 = 0.0226lo.ooo6 obtained 
by analyzing the cosmic microwave background radiation anisotropy. However, in high-redshift clouds, under 
conditions of low metallicity and low dust content, hydrogen may be incompletely molecularized even in the case 
of self-shielding. In this situation, the HD/2H2 ratio may not correspond to the actual D/H isotopic ratio. We 
have estimated the cloud molecularization dynamics and the influence of cosmological evolutionary effects on it. 

Key words, cosmology, primordial composition of matter, molecular clouds in early Universe, 
qso: J0812+3208, Q1331+170. 



1. Introduction 

Being the main coolant of the primordial gas, molecu- 
lar hydrogen (both H2 and HD) plays a central role in 
the creation of gas condensations and the formation of 
the first stars in the post-recombination Universe (see, 
e.g., Puy et al. 1993; Palla and Galli 1995; Lepp et al. 
2002; McGreer and Bryan 2008; Bromm et al. 2009). 
Observations of these molecules in high-redshift absorp- 
tion clouds allow the physical conditions existed in the 
early Universe to be determined. 

The fact that measuring the HD column den- 
sity makes it possible to independently determine 
the primordial D/H isotopic ratio is important too 
(Noterdaeme et al. 2008a; Ivanchik et al. 2010). Among 
the light nuclides produced by primordial nucleosynthe- 
sis 2 D, 3 T, 3 He, 4 He, 6 Li, 7 Li, and 7 Be, the primordial 
deuterium abundance is the most sensitive indicator of 
one of the key cosmological parameters - the baryon-to- 
photon ratio, r\ = n(,/fr 7 , or the corresponding baryon 
density of the Universe tlb (see, e.g., Sarkar 1996; Olive 
et al. 2000; Coc et al. 2005; Fields and Olive 2006). 

Until recently, the D/H ratio has been determined 
mainly only from H I and D I atomic lines in the absorp- 
tion spectra of quasars. However, such measurements 
run into a number of difficulties. The optical D 1 and 
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Hi spectra are almost identical; only the wavelengths 
of their lines are shifted by 0.027%. At the same time, 
the number densities of these atoms differ by four or 
five orders of magnitude. Therefore, if the Hi column 
density is low, then the D 1 lines are not seen at all. If, 
alternatively, the hydrogen column density is too high, 
then the H 1 lines are saturated, broadened, and overlap 
the Di lines (blend them). Moreover, the lines identi- 
fied as D 1 ones can in principle be produced by a small 
H I cloud that moves relative to the cloud being studied 
with a velocity of ~80 km/s, especially since there are 
actually many such clouds moving with different veloc- 
ities (the socalled Lyman-ev forest) on the line of sight. 
These factors may also be responsible for the significant 
spread in the D/H values obtained by this method. 

No difficulty with the identification of lines arises if 
the relative abundance of not the D 1 and H 1 atoms but 
the HD and H2 molecules is measured, because their 
spectra differ significantly and most narrow absorption 
lines do not overlap. However, the detection of molecules 
at high redshifts is a fairly rare event. For example, 
H2 molecules are observed only in ^10% of the DLA 
systems (Noterdaeme et al. 2008b) (damped Lyman-ev 
systems - those with a high atomic hydrogen column 
density, N H i > 10 20 cm- 2 , Wolfe et al. 2005) and only 
19 high-redshift H 2 systems have been detected to date. 
Only in three of them have HD lines been detected as 
well (the first system is Q 1232+0815 (Varshalovich et 
al. 2001; Ivanchik et al. 2010); the second system is 
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J 1439+1117 (Noterdaeme et al. 2008a); and the third 
system is J 2123-0050 (Malec et al. 2010)). The H 2 sys- 
tems are difficult to detect for two reasons. First, the 
molecular clouds are compact objects and, hence, the 
probability of their falling on the quasar-observer line 
of sight is low. Second, the H 2 absorption lines are 
very narrow and to detect them high-resolution spectra 
with high signal-to-noise ratio are needed. The chance 
of obtaining spectra with such a quality increased when 
large optical telescopes were put into operation. The HD 
molecular systems are even more difficult to observe, be- 
cause the number density of HD is several orders of mag- 
nitude lower than that of H2 . For the HD lines to be seen 
in a molecular cloud, the H 2 column density must be 
high enough, logNfj 2 > 19. This inequality may be con- 
sidered as a search criterion for HD systems. Therefore, 
out of the 16 remaining H 2 absorption systems, we gave 
attention to two systems with logNn 2 > 19. These are 
the absorption systems at z a b s = 2.626 in the spectrum 
of the quasar J 0812+3208 (logN H2 = 19.88, Jorgenson 
et al. 2009) and z a b s = 1.777 in the spectrum of the 
quasar Q 1331+170 (logN H2 = 19.65, Cu iet al. 2005). 
Having analyzed these systems, we detected two new 
HD systems (each of them consists of two subsystems) . 

2. Observational Data 

J 0812+3208. A spectrum of J 0812+3208 was taken 
with the HIRES echelle spectrograph at the Keck I tele- 
scope. This quasar was observed in the wavelength range 
3300-6200 A in 2005 under Ellison's and Wolfe's pro- 
grams (four 1-h and two 1.5-h exposures) and in 2007 
and 2008 under Prochaska's program (four 1-h and four 
40-min exposures) . The data are in free access from the 
Keck/HIRES arhciv^J We reduced and added these ex- 
posures using the MAKEE software packageQ specially 
developed by Barlow for Keck/HIRES. The total expo- 
sure time was 13.7 h, which allowed a signal-to-noise 
ratio of ~45 to be achieved at a resolution of ^46 000 
for the part of the spectrum containing H 2 and HD ab- 
sorption lines. 

Q 1331+170. A high-resolution spectrum of 
Q 1331+170 in the wavelength range 2280-3370 A, 
within which the molecular hydrogen lines at z = 1.777 
fall, was taken with the STIS spectrograph on the 
Hubble Space Telescope (HST) (program 7271 in 
1999 and program 9172 in 2002-2003, the principal 
investigator is Bechtold). The data are in free access 
from the HST archive^. The exposures were reduced 
and added using the CALSTIS procedures of the IRAF 
software package. The total exposure time was 14.5 h 
at a resolution of ~25 000 and the best signal-to-noise 
ratio of ~7. 

3. Data Analysis 

J 0812+3208. A DLA system was identified in the 
spectrum of J 0812+3208 (z cm = 2.701) in 2000 (White 
et al. 2000). The heavy-element abundance in this sys- 
tem was analyzed by Prochaska et al. (2003). Analysis 

1 http:/ /www2. keck, hawaii.edu/koa/koa. html 

2 http:/ /spider. ipac.caltech.edu/staff/tab/makee/index. html 

3 |http: / / archive.eso.org/archive/hst71 



of the C I absorption lines in this system (Jorgenson et 
al. 2009) revealed a cold cloud at redshift z a b s = 2.626. 
The upper limit on the cloud temperature T< 78 K was 
estimated from the Doppler parameter of neutral car- 
bon lines. Molecular hydrogen lines associated with the 
absorption system were also detected in the spectrum of 
this quasar. 

Q 1331+170. A DLA system with logN H2 =21.18 
was identified in the spectrum of Q 1331+170 (Carswell 
et al. 1975). The redshift of this system, z = 1.777, does 
not allow the part of the quasar spectrum with molecu- 
lar hydrogen lines to be obtained by ground based tele- 
scopes. Only HST observations allowed molecular hy- 
drogen lines to be identified in this system (Cui et al. 
2005). 

3.1. f/ 2 Column Densities 

J 0812+3208. We performed an independent analy- 
sis of the molecular hydrogen absorption system in the 
spectrum of J 0812+3208. A synthetic spectrum was 
constructed to determine the H 2 column density. The 
molecular hydrogen lines show the presence of two sub- 
systems in the spectrum, za = 2.626443(2) (subsystem 
A) and zb = 2.626276(2) (subsystem B), with a relative 
shift of ~14 km/s. This structure is most pronounced 
in the H 2 lines associated with the excited J = 2, 3, 
4 rotational levels of the X 1 !!^, v — ground state. 
The lines associated with the ground states of para- 
and orthohydrogen (J = and 1) are strongly satu- 
rated; it is hard to resolve the two subsystems in them 
and to determine the Doppler parameter b. Therefore, 
when determining the H 2 column density, we proceeded 
as follows. First, we analyzed the J = 2, 3, 4 levels. 
A local continuum was constructed for each of the se- 
lected lines (a total of ^50 lines belonging to the Lyman 
(from LO-0 to L17-0) and Werner (from WO-0 to W4-0) 
bands were used). The Doppler parameter was assumed 
to be the same for all levels. It should be noted that the 
parameter b may grow with increasing rotational level 
number (see, e.g., Balashev et al. (2009) and references 
therein). Here, we disregarded this possibility, because 
our main objective was to determine the total hydro- 
gen column density, which for this system is determined 
mainly by the ground rotational levels whose lines are 
strongly saturated and, as a result, are insensitive to 
the parameter b. Once the values of the redshifts and 
Doppler parameters for the two subsystems had been 
obtained, we fixed these values and constructed a syn- 
thetic spectrum for the ground H 2 levels. The synthetic 
H 2 spectrum is shown in Fig. [1] and the parameters ob- 
tained from our analysis are presented in Table [TJ For 
subsystem A, our H 2 results are in good agreement with 
those of Jorgenson et al. (2009). However, the latter au- 
thors determined the Doppler parameter bn 2 from the 
Doppler parameter of neutral carbon bci, while we used 
it as an independent parameter being fitted. 

For subsystem B, the derived column densities of H 2 
molecules in the J = and 1 ground states have large 
errors, because the column density in subsystem B is 
an order of magnitude lower than that in subsystem A. 
Therefore, the saturated lines are centered near za due 
to the small relative shift between subsystems A and B, 
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Figure 1. Synthetic H2 spectrum of the absorption system at z = 2.626 fitted into the observed spectrum of 
J 0812+3208 (HIRES/Keck). 



while subsystem B makes only a slight contribution to 
the asymmetry of the Lorentz wings. 

Fig.[2]shows the relative populations of H2 rotational 
states in the two subsystems. Using the ratio of the 
ortho-H2 and para-H2 (J = and 1) column densities, 
we estimated the gas temperature in both subsystems: 



T#i = 48±2K and 



50±?«K. 



Q 1331+170. In contrast to the analysis of the H2 
absorption system in Cui et al. (2005), where the au- 
thors fitted only one component (z = 1.776553(3)) into 
the spectrum, we constructed a synthetic spectrum by 
assuming the absorption system to consist of two com- 
ponents (z c = 1.77637(2) and z D = 1.77670(3)). Two 
components with a relative shift of ^35 km/s are clearly 
seen in the H2 lines originating from the J = 4 rotational 
level (Fig. [3]) and in the HD lines we identified (see 
Fig. [5] below). The redshifts of both subsystems were 
determined by analyzing the H2 (J = 4) lines. Since the 
lines associated with the J = and 1 levels are strongly 
saturated for both subsystems, the column density of 
molecules at these levels can be obtained with a good ac- 
curacy. This cannot be said for the lines associated with 
the 3 = 2 and 3 levels, which fall into the logarithmic 
part of the curve of growth. As a result, the total molec- 
ular hydrogen column density for the two subsystems is 
logNg 2 = 19.43 + 0.10 andlogNg 2 = 19.39 + 0.11. 

3.2. HD Column Densities 

J 0812+3208. To determine the HD column density, 
the Doppler parameter, and the redshift, we constructed 
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Figure 2. Populations Nj of H2 levels (normalized to 
their statistical weights gj) versus level energy Ej. The 
red and blue circles correspond to subsystems A and 
B in the spectrum of J 0812+3208. The straight lines 
corresponds to Tqi, the excitation temperatures of the 
J = and J = 1 levels. 



a synthetic spectrum for each of the two subsystems (A 
and B). The parts of the quasar spectrum normalized to 
the local continuum and the synthetic spectrum for HD 
lines are shown in Fig. @] In subsystem A, we identified 
the HD lines corresponding to the transitions not only 
from the J — ground state but also from the J = 1 
rotational level. This allowed us to estimate the number 
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Table 1. Results of our analysis of the absorption system at z = 2.626 in the spectrum of J 0812 + 3208 



Subsystem A 
z A = 2.626443(2) 



Subsystem B 
z B = 2.626276(2) 







log(N) 


b, km/s 






log(N) 


b, km/s 


H 2 


J=0 


19.83 ±0.05 


0.81 ±0.10 


H 2 


J=0 


18.71 ±0.45 


1.55 ± 0.11 




J=l 


19.25 ±0.02 


0.81 ±0.10 




J=l 


18.19 ±0.19 


1.55 ± 0.11 




J=2 


16.47 ±0.10 


0.81 ±0.10 




J=2 


16.79 ±0.10 


1.55 ± 0.11 




J=3 


15.15 ±0.25 


0.81 ±0.10 




J=3 


14.62 ± 0.09 


1.55 ± 0.11 




J=4 


13.95 ±0.12 


0.81 ±0.10 




J=4 


13.39 ± 0.14 


1.55 ± 0.11 




Ej 


19.93 ± 0.04 
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18.82 ± 0.37 
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J=0 


15.70 ±0.07 


0.7 ±0.04 
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J=0 


12.98 ± 0.22 


< 3.6 1 




j=i 


13.77 ±0.15 


0.7 ±0.04 
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15.71 ±0.07 
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1 - The observed lines are insensitive to the parameter b, because they are located at the linear part of the 
curve of growth. 
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Figure 3. Synthetic H2 spectrum of the absorption system at z = 1.777 fitted into the observed spectrum of 
Q 1331±170 (STIS/HST). The two components of the synthetic spectrum are indicated by the dotted and dashed 
lines. These components are well resolvable for the lines originating from the J = 4 rotational level (bottom). 



density in the cloud, n = 541^2 cm~ 3 (for more detail, 
see below). As far as we know, this is the first reliable 
detection of HD lines associated with the J = 1 level. 
This proved to be possible due to the high HD column 
density and good spectrum quality. 

The HD molecular lines of the Lyman and Werner 
series are present in the spectrum up to L18-0 and W4-0 
inclusively. In general, the HD wavelengths and oscilla- 
tor strengths for Lz/' — and Wv" — lines with large v" 
are not known well enough. The laboratory wavelengths 
were measured with a high accuracy, 8X/X ~ 5 x 10 , 
only for the bands up to L9-0 and WO-0 (Ivanov et 
al. 2008). For shorter wavelength transitions, we used 
the laboratory data (Ivanov et al. 2010) kindly pro- 
vided to us by Prof. W. Ubachs. The oscillator strengths 
were taken from the calculations by Abgrall and Roueff 
(2006); the wavelengths are also given there. However, 
comparison of the theoretical (Abgrall and Roueff 2006) 



and laboratory (Ivanov et al. 2008) values showed that 
the calculated wavelengths for some lines differ signifi- 
cantly from the laboratory ones. For example, the cal- 
culated wavelength of the L11-0R(0) line differs from 
the laboratory value by ^8 km/s and does not fit into 
the general analysis. Using the data from Ivanov et al. 
(2010), we were able to determine the HD column den- 
sity with a good accuracy. 

The results of our determination of the parame- 
ters for the identified HD molecular lines in the spec- 
trum of J 0812±3208 are presented in Table □ The red- 
shifts of the HD lines in subsystems A and B coincide, 
within the error limits, with those determined by an- 
alyzing H2 lines. To determine the HD (J = 0) col- 
umn density in subsystem A, we chose the R(0) lines 
of the following bands: L0-0, L3-0, L4-0, L5-0, L8-0, 
Lll-0, L18-0, W0-0, W4-0. The 1, 2, and 3cr confi- 
dence regions for Nhd and bnD are shown in Fig. [SJ 
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Figure 4. Synthetic HD molecular spectrum of the absorption system at z = 2.626 in the observed spectrum of 
J 0812+3208. Two components are seen for the transitions from the J — ground state. The R(l), P(l), and Q(l) 
lines of the transitions from the excited J = 1 rotational level are also shown. 



We determined the HD column density in subsystem 
A, log Njjd = 15.71 ± 0.07. This is the highest value 
of Nhd measured in the quasar spectra. The derived 
Doppler parameter bfjD = 0.70 ± 0.04 km/s agrees 
(within the error limits) with bci = 0.33 ± 0.05 km/s 
and bn 2 = 0.81 ± 0.10 km/s (if the thermal broadening 
is dominant, then the Doppler parameters are related as 
the square root of the ratio of the element masses). 



Number density in the subsystem 
J0812+3208A. The relative populations of the 
HD J = and J = 1 rotational levels determined 
in the absorption subsystem J0812+3208A allow the 
number density of the molecular cloud to be estimated. 
This requires considering the level population balance 
equations. The J = 1 level can be populated by 
the direct radiative transition from the J = level, 
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Figure 5. Confidence regions of the parameters Nhd 
and t>HD for the J = level of subsystem A in the spec- 
trum of J 0812+3208. The minimum corresponds to the 
column density log Nhd = 15.70 and the Doppler pa- 
rameter t>HD = 0.7 km/s. 



radiative pumping through excited electronic states, 
and collisions. Typically, the particle number density in 
a diffuse cloud is n^ 10 — 500 cm~ 3 . At such number 
densities and the Galactic mean intensity of ultraviolet 
radiation, collisions will be the dominant J = 1 level 
population channel for the bulk of the cloud. Indeed, 
at such a high measured column density in the cloud, 
log Nhd = 15.7, the rate of the Q—tl population by ra- 
diative pumping inside the cloud falls by more than two 
orders of magnitude through selfshielding. The direct 
radiative 0— >1 population can also be neglected. The 
dominant channel of the 1— »0 transition is spontaneous 
relaxation characterized by the radiative transition 
probability Aio- Thus, we can estimate the number 
density from the formula 

N J= i A 01 
Nj= C01 

The main collisional parter of HD for a molecular 
cloud will be H2; the values of the collisional prob- 
ability coefficient C01 were taken from Flower et al. 
(2000). Using a temperature estimate in the cloud of 
T£ = 48 ± 2 K, we obtained the following estimate 
for the number density in the cloud: n= 541^2 cm -3 it 
agrees with typical values for diffuse clouds (Rachford et 
al. 2002). Note, however, that if the cloud is not molcc- 
ularized, then the main collisional partner is atomic hy- 
drogen and the number density estimate doubles. 

Q1331+170. The HD molecular lines are present in 
the spectrum up to the L18-0 Lyman and W4-0 Werner 
bands, respectively (some of these lines are shown in 
Fig. [5]). Because of the low signal-to- noise ratio, we were 
able to use only the L3-0R(0), L4-0R(0), L7-0R(0), and 
L8-0R(0) lines to determine the HD parameters. Two 
subsystems whose redshifts coincided, within the error 
limits, with those determined by analyzing H2 are seen 
in the chosen lines. The HD column densities determined 
by synthetic spectrum construction were found to be 
logNgpj = 14.83 ± 0.15 and logNg D = 14.61 ± 0.20. 



3.3. HD/2H 2 Ratio 

Until now, HD lines (at high redshifts) have been iden- 
tified only in three of the 19 H2 absorption systems. 
Figure [7] shows the HD, H2 and Dl, Hi column densi- 
ties measured in interstellar clouds of our Galaxy and 
absorption systems of quasars. Determining the column 
densities of HD and H2 or D 1 and H 1 at high redshifts 
allows the primordial D/H isotopic ratio and, hence, Qb 
to be estimated. The oblique straight lines in Fig.[7]indi- 
cate the column densities that correspond to the mean 
D/H and HD/2H2 ratios for the Galactic and extra- 
galactic systems. In our Galaxy, the D/H ratio mea- 
sured from the column densities of atomic D I and H 1 
lines (Linsky et al. 2006) is systematically lower than 
its mean value obtained by analyzing the quasar spectra 
(Pettini et al. 2008). In principle, this can be explained 
by the astration of deuterium in stars. We also see that 
there is a significant difference between the values of 
D/H and HD/2H2 measured in our Galaxy. This may be 
because in contrast to H2, the HD molecules are not al- 
ways shielded from ultraviolet radiation and, hence, deu- 
terium is molecularized to a lesser extent. Another ex- 
planation for the shortage of the HD molecular fraction 
can be complex chemistry of molecular clouds, where 
deuterium can effectively enter into other, more com- 
plex molecules, H 2 0, NH 3 , HCN, polyaromatic hydro- 
carbons, etc. No significant difference is observed for the 
D /H ratios in the spectra of quasars measured from the 
molecular and atomic components (Fig.[7]). However, the 
statistics for quasar spectra is much poorer than that for 
systems in our Galaxy. 

J 0812+3208. The ratios of the HD and H 2 column 
densities in the spectrum of J 0812+3208 were found 
to be HD/2H 2 = 2.97±g;|g x 10~ 5 and HD/2H 2 = 

7.08t4° 2 g 5 x 10~ 7 for the absorption systems A and B, 
respectively. The HD/2H2 ratio in subsystem B is defi- 
nitely unsuitable for the D/H estimate, because the HD 
column density in this subsystem (log Nhd = 12.98) is 
two orders of magnitude lower than the value necessary 
for the self-shielding of molecules and, as a result, deu- 
terium may be incompletely molecularized. According 
to the calculations by Le Petit et al. (2002), if HD 
and H2 are self-shielded from ultraviolet radiation, then 
both components arc completely molecularized and the 
universal ratio D/H = HD/2H2 is established in the 
cloud. In our case, only subsystem A can be completely 
shielded in HD, because it has a relatively high HD col- 
umn density: log Nhd = 15.7. Besides, this subsystem 
has a low metallicity (relative to the solar one), < — 1.0 
(Prochaska et al. 2003). This means a low level of deu- 
terium astration in stars and makes it possible to es- 
timate the primordial D/H ratio. On the other hand, 
a low metallicity also means a low dust content; since 
molecular hydrogen is formed mainly on dust, this leads 
to a lower formation rate of H2 and HD and, in con- 
trast to the calculations by Le Petit et al. (2002), the 
cloud molecularization may be incomplete. Since there is 
an additional, fairly efficient formation channel for HD, 
H2+D+ — > HD + H + , at a low H2 formation rate on dust 
this can lead to a significant increase in the HD/2H2 
ratio and it will not correspond to the true D/H iso- 
topic ratio (see, e.g., Ferlet et al. 2000). Besides, at a 
low dust content, the process of gas molecularization in 
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Figure 6. Synthetic HD spectrum of the absorption system at z = 1.777 fitted into the observed spectrum of 
Q 1331+170. The R(0) lines of the transitions from the J = ground state, in which two components are clearly 
seen, are shown. 



the cloud slows down. Therefore, the matter in the ob- 
served clouds may not have time to reach a steady state 
in chemical equilibrium and the HD/2H2 ratio can then 
be different for clouds of different ages (see the Section 0] 
and Fig.©. 

Q 1331+170. The ratios of the HD and H 2 column 
densities in the spectrum of Q 1331+170 were found 
to be HD/2H 2 = 1.24±g;|o x 10~ 5 and HD/2H 2 = 
0.831 

0.33 x 10 5 for the absorption systems C and D, 
respectively. Note that the large errors in the parame- 
ters compared to the previous system are attributable 
to a poorer spectrum quality. 



3.4. Peculiarities of the HD/H2 Absorption Systems 

J 2123-0050. According to Malec et al. (2010), the ab- 
sorption system in the spectrum of J 2123-0050 has the 
highest HD/2H2 ratio, ~ 7 x 10~ 5 . However, in our opin- 
ion, this system requires an additional study of the col- 
umn densities and its physical parameters, because the 
ortho-to-parahydrogen ratio obtained by these authors, 
~ 16:1, exceeds considerably the ratio 3:1, appears very 
strange, and has never been observed before. 

J 1439+1117. The metallicity of this system is 
nearly solar and HD/2H2 = 1.5 x 10~ 5 is close to the 
D/H ratio measured in our Galaxy (Noterdaeme et al. 
2008a). 

Q 1232+0815, J 0812+3208. This system has a 
low metallicity, ~ —1.5 (Balashev et al. 2010; Prochaska 
et al. 2003), and high HD column densities. This makes 
the derived values for these HD/2H2 systems a good 
estimate of the primordial D/H ratio. 

Q1331+170. The metallicity of this system is also 
low, ~ —1.5 (Prochaska and Wolfe 1999). However, in- 
complete shielding in HD is quite possible for the two 
derived HD/H2 systems, because the HD column density 
is low. Besides, the spectrum has a low signal-to-noise 
ratio and, hence, the estimate of the primordial D/H 
ratio based on these systems is ambiguous. 

The results for all the systems are summarized in 
Table H 



4. Cloud Molecularization Model 

To estimate the influence of incomplete moleculariza- 
tion on the HD/H2 ratio, we calculated the chemical 
evolution of an HD/H2 molecular cloud (this model 
will be presented in detail in Balashev et al. (2010)). 
Figure [8] presents the results of our calculations of the 
chemical evolution of a typical interstellar cloud with 
a temperature T = 80 K and a total hydrogen num- 
ber density n = 100 cm~ 3 . We consider a cloud com- 
posed of H, He, D, their ions, electrons, H2 and HD 
molecules, and some amount of dust on which these 
molecules are efficiently formed. The chemical reaction 
rates were taken from Standi et al. (1998) and Glover 
and Jappsen (2007). We also took into account more 
complex molecules, H +3 , H2D+, HcH + , HcD + , etc., but 
their influence on the HD/H2 evolution of interest to us 
proved to be insignificant in the case of a diffuse molec- 
ular cloud. The cloud was assumed to be in the field of 
external ultraviolet radiation, which leads to an efficient 
dissociation of molecules in the cloud surface layers. Its 
intensity was taken to be equal to the mean intensity 
of ultraviolet radiation in our Galaxy. The function of 
cloud self-shielding from ultraviolet radiation was taken 
from Draine and Bertoldi (1996). We also assumed that 
the cloud was irradiated by cosmic rays, which leads 
to gas ionization that provides ionmolecular reactions. 
The cosmic-ray ionization rate was taken to be equal 
to its Galactic mean value. The dust content was taken 
to be 1 (Fig. 8, left) or 0.01 (Fig. 8 right) of the mean 
dust content in our Galaxy. These values were chosen 
so as to qualitatively explain the difference between the 
observed clouds in our Galaxy and high-redshift DLA 
systems. The total hydrogen column density in the sys- 
tem Ntot = Nhi + 2Nh 2 was chosen from the condi- 
tion Nh 2 — 10 20 cm -2 , which corresponds to the column 
density for the absorption subsystem A of the quasar 
J 0812+3208. 

We see that this model can explain the observational 
data, provided that the cloud age is greater than 10 7 yr. 
At the Galactic dust content, the cloud reaches almost 
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log 2N(H 2 ) 

Figure 7. D I, Hi and HD, H2 column density measurements. The circles indicate the results of HD and H2 
measurements in the quasar spectra (the red circles correspond to the data obtained in this paper); the squares 
represent the results of HD and H2 measurements in our Galaxy (Snow et al. 2008; Lacour et al. 2005); the blue 
and light blue squares indicate the points for which the HD column density is higher and lower than logNnD = 
15, respectively. The solid straight line indicates the mean HD/2H2 ratio determined from points with logNnD > 
15; the green stripe indicates the results of Hi and Di measurements taken from Linsky et al. (2006) (the mean 
D/H ratio for these values is indicated by the dotted line). The D/H ratio in the quasar spectra measured from 
H I and D I lines (Pettini et al. 2008) is indicated by the yellow stripe. 
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Figure 8. Time variation of the HD/2H2 abundance ratio (top panels) and the molecular fraction (bottom panels) 
during the evolution of an interstellar cloud with a number density n= 100cm -3 and T = 80 K. The calculations 
are presented for the Galactic dust content (left panels) and 0.01 of the Galactic dust content (right panels). The 
vertical line marks the characteristic dynamical lifetime of clouds, ~3x 10 7 yr. 
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QSO 



Njt cm 



N N d, cm 2 



Instrument 



References 1 



Q 1232+0815 
J1439+1117 2 
J 2123+0050 
Q 0812+3208 
Q 0812+3208 
Q 1331+170 
Q 1331+170 



2.33714(3) 
2.41837 
2.0593276(5) 
2.626443(2) 
2.626376(2) 
1.77637(2) 
1.77670(2) 



4.78±g;|| x 10 19 



2.40±°;*j x 10 19 



3-68±o.3? x 10 17 



8.54± ;?| x 10 19 
6.74±|;|1 x 10 18 
2.66±°;|1 x 10 19 



2.45l°;g x 10 19 



3.39±o.l x 10 15 



7.18l°;« x 10 1 



5.89lg;S x 10 1 



5.07±g;?l x 10 1 



4.07±IH x 1Ql4 



UVES 
UVES 
KECK 
KECK 
KECK 
HST 
HST 



1, 2 
3 
4 
5 
5 
5 
5 



1 Refs.: 1 Varshalovich et al. (2001), 2 Ivanchik et al. (2010), 3 Noterdaeme et al. (2008a), 4 
Malec et al. (2010), 5 this paper. 

2 In this system, three components are seen in the HD lines; however, since the H2 column densities 
cannot be determined for these components, the authors provide only an integrated value. 



2N H2 



1 



complete molecularization in H 2 , ih 2 = 2 n„ +n hi 
(Figs. 8 bottom left); the HD molecular fraction {iiiD = 
nJd+No f ) 1S l° wer ' because HD is shielded more weakly 
against ultraviolet radiation than H 2 . The situation is 
different at a low dust content. The hydrogen molecular 
fraction is low (fn 2 ,(hd < 0.05); the ion reaction of HD 
formation becomes more efficient than the formation of 
HD on dust and, as a result, the molecular fraction of 
HD can be higher than that of H 2 (Figs. 8 bottom right). 

It is important to note that HD/H 2 depends signif- 
icantly on the cloud age at a low dust content (Figs. 
8 top panels). In addition, the HD/H 2 ratio also de- 
pends on other factors, namely, the ultraviolet radia- 
tion background, the cosmic-ray intensity, the gas num- 
ber density, and the total column density. Therefore, 
the HD/2H 2 ratio may not correspond to the actual 
D/H isotopic ratio and its estimation requires devel- 
oping a model for the cloud molecularization dynamics 
(Balashev et al. 2010). Note that for the two most satu- 
rated (in HD) systems (J 0812+3208A and Q 1232+082, 
see Fig. 7), the HD/2H 2 ratios turned out to be coinci- 
dent (within the error limits) and close to the primor- 
dial D/H determined by analyzing the cosmic microwave 
background radiation anisotropy (Komatsu et al. 2010). 
It is hoped that increasing the statistics of high-redshift 
molecular clouds will clarify this situation. 

In molecular cloud models, the D/HD transition is 
assumed to occur at larger penetration depths of the 
emission into the cloud than the H/H2 transition (see, 
e.g., Le Petit et al. 2002). However, under conditions of 
low dust content, a situation where the D/HD transition 
in the cloud will begin more early than the H/H 2 tran- 
sition is possible. This will occur via the ion reaction 
of HD formation under incomplete H 2 molecularization 
conditions. 



At a low dust content and high column densities, 
Nh 2 > 10 21 cm~ 2 , on long time scales (> 10 8 yr), the 
cloud reaches complete molecularization in H 2 and HD 
and the HD/2H 2 ratio will corresponds to the D/H ratio. 
However, there will exist a long time interval (~ 10 7 yr) 
in the cloud evolution in which the HD/2H 2 ratio will 
exceed considerably the D/H ratio. 



5. Conclusions 

We identified HD molecular lines in the absorption sys- 
tems at z a b s = 2.626 and z a b s = 1.777 in the spectra of 
the quasars J 0812+3208 and Q 1331+170, respectively. 
Each system consists of two clearly resolvable compo- 
nents. The HD and H 2 column densities were deter- 
mined for all four subsystems (Table 2). 

Since only one of the four identified subsystems, 
at z = 2.626443(2), in the spectrum of J 0812+3208 
has an HD column density exceeding its critical value 
for self-shielding, one would expect the HD/2H 2 ra- 
tio in the steady-state model of a completely molecu- 
larized cloud to characterize the primordial D/H iso- 
topic ratio. Under these conditions, D/H= (2.97lg;gg) x 
10~ 5 , which corresponds to the baryonic matter den- 
sity n h h 2 = 0.0205^ ° 2 ° 25 . This value is in good agree- 
ment with fibft 2 = 0.0226Jq ° ° 6 06 obtained by analyzing 
the cosmic microwave background radiation anisotropy 
(Komatsuet al. 2010). 

However, the cloud molecularization dynamics and 
the influence of cosmological evolutionary factors re- 
lated to a change in dust content, ionizing radiation 
background, and cosmic-ray intensity on it can lead to 
a systematic shift in the estimate of the primordial D/H 
isotopic ratio. 

We made preliminary estimates of the effects related 
to the cloud molecularization dynamics in the early 
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Universe. Their allowance was shown to be necessary 
to properly determine the primordial D/H ratio. 

Increasing the statistics of HD/H2 molecular sys- 
tems will allow progress to be made both in understand- 
ing the chemical and isotopic evolution of molecular 
clouds as well as for determining the primordial D/H 
isotopic ratio. Significant progress in solving this prob- 
lem may be expected when future extremely telescopes 
will be put into operation. 
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Figure 9. Synthetic H2 spectrum of the absorption system at z = 1.777 fitted into the observed spectrum of 
Q 1331+170 (STIS/HST). 



